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ABSTRACT 

We have obtained multi-epoch, high-resolution spectroscopy of 218 candidate low- 
mass stars and brown dwarfs in the young clusters around <j Ori and A Ori. We find 
that 196 targets are cluster members based on their radial velocity, the equivalent 
width of their NaI8200 lines and the spectral type from their TiO band strength. We 
have identified 1 1 new binary stars among the cluster members based on their variable 
radial velocity and an additional binary from the variation in its line width and shape. 
Of these, 6 are double-lined spectroscopic binaries (SB2) where the components of the 
binary are of comparable brightness. The others are single-lined binaries (SBl) in 
which the companion is faint or the spectra of the stars are blended. There are 3 
narrow-lined SBl binaries in our sample for which the companion is more than 2.5 
magnitudes fainter than the primary. This suggests that the mass ratio distribution 
for the spectroscopic binaries in our sample is broad but that there may be a peak 
in the distribution near q = l. The sample covers the magnitude range Ic = 14- 18.9 
(mass « 0.55 — O.OSMq), but all of the binary stars are brighter than Ic = 16.6 
(mass K, 0.12Mq) and 10 are brighter than Ic = 15.5 (mass ~ 0.23Mq). There is a 
significant lack of spectroscopic binaries in our sample at faint magnitudes even when 
we account for the decrease in sensitivity with increasing magnitude. We can reject 
the hypothesis that the fraction of spectroscopic binaries is a uniform function of Ic 
magnitude with more than 99 percent confidence. The spectroscopic binary fraction 
for stars more massive than about O.IMq (Ic < 16.9) is /bright = 0-095lQ;Q2g. The 90 
percent confidence upper limit to the spectroscopic binary fraction for very low mass 
(VLM) stars (mass < O.IMq) and brown dwarfs (BDs) is /faint < 7.5 percent. The 
hypothesis that /bright and /faint are equal can be rejected with 90 percent confidence. 
The average detection probability for our survey is 50 percent or more for binaries 
with separations up to 0.28 au for stars with Ic < 16.9 and 0.033 au for the fainter 
stars in our sample. We conclude that we have found strong evidence for a change 
in the fraction of spectroscopic binaries among young VLM stars and brown dwarfs 
when compared to more massive stars in the same star-forming region. This implies a 
difference in the total binary fraction between VLM stars and BDs compared to more 
massive stars or a difference in the distribution of semi-major axes, or both. 

Key words: binaries: spectroscopic - stars: low-mass, brown dwarfs. 



1 INTRODUCTION 



The origins of very low-mass stars (VLMS) and brown 
dwarfs (BD) are proving difficult to understand, despite be- 
ing more common than stars of higher mass. Ideas include 
ejection from protostellar aggregates (|Reipurth fc Clarke! 
[20011 fo rmation within convergent flows generated by tur- 



bulence (|Padoan fc Nordlund 2004), the photo- erosion of 
pre-stellar cores i Whitworth fc Zinneckeil |2004 ) or frag- 



mentation within the outer pa rts of circumstellar discs 
jWhitworth fc Stamate"iio3l2006l ). 



The frequency and separation distribution of binary sys- 
tems is an important constraint on the likely formation pro- 
cess. There is strong evidence that the binary properties 
of the lowest mass stars and brown dwarfs are quite dif- 
ferent to those of hi gher mass objects (see the review of 
iBurgasser et al.ll2007l ). Resolved imaging of nearby VLMS 
and BDs in the field show that about 15-20 percent of 
systems are binaries with separations greater than 1-2 au, 
but that very few have separa t ions g reater than 20 au (e.g. 
IClose et all l2003l : IBouv et al.l |2003| ). This contrasts with 
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higher mass stars where overaU binary frequencies are 30-60 

J iercent, with a much broader spread of possible sep arations 
Fischer fc Marcvlll99i iDuauennov fc Mavod ll99ll) . 

A missing part of the picture is how many VLMS 
and BDs are binary systems with separations less than 
about 1 au, where the im aging observations cannot reach. 
iMaxted fc Jeffries! l|2005l ) used previously published radial 
velocity (RV) results to show that an overall binary fre- 
quency (at all separations) of 32-45 percent was needed 
to explain the presence of several binaries detected by RV 
variations, with most of them at small separations. This 
high frequency found support from a small RV survey by 
iKurosawa et al.l l|2006l ) who found a frequency of 11 -40 per- 
cent for separations less than 0.1 au i n a yo ung association. 
On the other hand iBasri fc ReinersI l|2006l ) found few RV 
variables in their survey of field VLMS/BDs and concluded 
that the overall binary frequency was 16-36 percent, with 
very few binaries at separations below 1 au 

The situation is unresolved and it is quite likely that 
the apparent discrepancies between these various authors 
arise from biases within the samples considered; differences 
in analysis technique and that in terms of absolute numbers, 
very few short-period VLMS/BD binary systems have yet 
been found, thus limiting the statistical precision possible. 

In this paper we present the results of an RV sur- 
vey of a large number of low-mass stars and brown 
dwarfs in the a Ori and A Ori clusters. These clusters 
are young and nearby (3-5 Myr, 330-4 50 pc) and contain 
large populations of VLMS and BDs Js^iar et all [iooi : 
iBarrado v Navascues et al.1 12004| : iKenvon et al.l 120051 ). We 
have observed more than 200 objects using the flames 
multi-fibre spectrograph on the VLT-Kueyen telescope to 
measure radial velocities at several epochs and searched for 
close binary systems at a range of masses. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

We have used the flames multi-object spectrograph 
jPasauini et al.l 120021 ) on ESO's VLT UT2 (Kueyen) tele- 
scope to obtain multi-epoch, high-resolution spectroscopy 
of 218 faint stars in the clusters around a Ori and A Ori. 
The instrument is capable of providing spectra of up to 130 
targets in one setting over a field of view 25 arcmin in di- 
ameter. 

We sele cted targets around a Ori from the photometric 
catalogue of iKenvon et al.l (|2005l ) . We included all stars in 
the magnitude range 14 < Ic < 19 with the correct Ic mag- 
nitude and (Rc — Ic) colour to be a cluster member in the 
input catalogue for the fibre allocation process, irrespective 
of any other membership information that may have been 
available. For A Ori all the candidate memb e rs ide ntified in 
the catalogue of iBarrado v Navascues et al.l (|2004l ) were in- 
cluded in the input catalogue for the fibre allocation process. 
Fibres were allocated in 6 fields, 2 fields near A Ori and 4 
fields near a Ori. The central position of each field and the 
number of targets for which useful spectra were obtained 
are given in Table [T] The position of the targets in the Ic 
v. (Rc — Ic) colour-magnitude diagram is shown in Fig. [1] 
Those fibres that could not be allocated to targets were used 



Table 1. Log of observations. The number of stars observed with 
the GIRAFFE spectrograph on each date for which useful spectra 
(signal-to-noise ratio > 5) were obtained is given in the column 
headed N . The date of observation is given in the fifth column. 
The standard deviation of the radial velocities of the night sky 
emission line spectra obtained at each epoch, Usi^y, arc given in 
the final column. 



Field 


Position 
Q!,72()0() <5j20()0 


N 


Date 
(UTC) 


''"sky 

(kms-1) 


A Oril 


05 35 58.5 


4-09 51 21 


34 


20051014 


0.66 








34 


20051117 


0.59 








35 


2005 12 30 


0.67 


A Ori 2 


05 34 44.1 


-f09 5144 


49 


20051014 


0.64 








49 


2005 11 27 


0.99 








50 


2005 11 28 


0.63 








50 


2005 12 30 


0.58 


a Oril 


05 40 28.0 


-0215 35 


26 


2005 11 13 


0.60 








22 


20051117 


0.60 








26 


2006 01 08 


0.64 








26 


2006 01 13 


0.92 


o- Ori 2 


05 4001.6 


-02 40 02 


17 


2005 11 13 


0.67 








29 


2005 12 07 


0.65 








28 


2006 01 13 


0.74 


(7 Oris 


05 38 26.9 


-02 4121 


56 


2005 11 13 


0.64 








55 


2005 11 24 


0.67 








56 


200512 07 


0.56 








55 


2006 01 13 


0.60 


a Ori 4 


05 38 23.8 


-0214 59 


33 


20051111 


0.59 








34 


20051111 


0.58 








34 


2005 12 07 


0.62 








32 


2006 01 13 


0.60 



to obtain simultaneous background spectra for sky subtrac- 
tion. The two fields around A Ori overlap so 10 stars were 
observed in 2 fields. For most stars we have obtained 3-4 
spectra with the exception of 10 stars observed in 2 fields 
for which we typically have 6-7 spectra. The baseline of 
the observations is at least 24 days for all stars and is more 
typically 60-70 days. The sensitivity of our survey drops 
rapidly for binaries with semi-major axes a ^ 0.25 au as a 
result of the baseline of our observations. 

Light from the fibres was fed to the giraffe spectro- 
graph operated in a high resolution mode with the H836.6 
echelle grating. A filter was used to select light from the 
6**^ echelle order which covers the wavelength region 8073 - 
8632 A. The resolving power of the spectra is R ~ 16000. 
Spectra were obtained in service mode by ESO staff on 
the dates given in Table [T] The exposure time in each case 
was 2750 s. The seeing during the exposures was typically 
0.92arcsec but varied from 0.48arcsec to 1.45arcsec. 

We also took advantage of the possibility to use up to 8 
fibres from fl ames to obtain spe ctra with the UVES echelle 
spectrograph (jPekker et al.|[2000l ) at the same time as the 
GIRAFFE observations. We used the CD#4 echelle grating 
to obtain spectra at a resolution of i? ~ 47000 covering the 
same spectral range as the giraffe spectra for one or two 
bright, early-type stars in the field. These spectra were used 
to calibrate the telluric absorption in the giraffe spectra. 
Fibres not allocated to bright stars were used to obtain si- 
multaneous spectra of the night sky. 
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Figure 1. Colour - magnitude diagrams for stars around tr Ori (left panel) and A Ori (right panel) with the stars we have observed 
highlighted (open circles). Other symbols in the left panel indicate the membership status for stars from lKenvon et al.l (2005) as follows: 
filled circles - me mber; triangles — uncertain member ship; crosses — non-member. This membership status was ignored during the fibre 
allocation process. iBarrado v Navascues et al.l | |2004) only provide data for stars they consider to be candidate members of the A Ori 
cluster so we are not able to show that position of other stars in the right panel. 
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Figure 2. Example of a synthetic telluric absorption spectrum from a model terrestrial atmosphere fit to the telluric absorption in a 
section of one of our UVES spectra. The residuals from the fit of the synthetic absorption spectrum plus a low order polynomial are shown 
below the UVES spectrum and the synthetic absorption spectrum. Note that the residuals include absorption lines from the star, e.g., the 
Nal 8200A doublet. 



2.2 Reduction of the spectra 

There is strong fringing in the images produced by the GI- 
RAFFE spectrograph when operated at the wavelengths we 
have used for our observations. This makes the extraction of 
the spectra problematic. For brighter stars it is possible to 
use the spectra and their associated errors extracted from 
the data automatically by the ESQ pipeline provided by the 
observatory. These spectra are extracted from the images by 
summing the pixels in each row within a given range around 
the central position of each fibre. The disadvantages of us- 



ing these spectra are that there is no cosmic-ray rejection 
included in the extraction process and that the signal-to- 
noise ratio drops rapidly as the spectra become fainter. 

For fainter stars we tried using the giraffe Base- 
Line Data Reduction Software (girBLDRS) version 1.13.1 
(Blecha et al. 2000) to perform optimal extraction of the 
spectra. This maximizes the signal-to-noise in the result- 
ing spectra by weighting the pixels according to th e vari- 
ance of each pixel and a model of the spatial profile l|Hornel 
Il986l ). Some sort of weighted extraction is essential to pro- 
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Figure 3. Example spectra of apparently single stars from our sample. The spectra of each star have been normalized by a low-order 
polynomial and oflFset for clarity and are labeled by name and I-band magnitude. 



duce usable spectra for the faintest stars in our dataset. 
It is also possible to identify pixels affected by cosmic-rays 
in the images and exchide them from the extracted spec- 
tra by comparing the spatial profile of the spectra with a 
model. Unfortunately, we were not able to achieve the same 
level of stability in the radial velocities for data extracted 
with girBLDRS as the normal pipeline extraction. We be- 
lieve that this is due to girBLDRS separately optimizing the 
profile used to extract the object spectra, the flat-field spec- 
tra and the thorium-argon arc spectra. The result is that 
the profiles used to extract the object, arc and flat spectra 
are slightly different. This optimization is used to account 
for small shifts in the positions of the spectra on the de- 
tector. The fringe-pattern can change the effective detector 
efficiency by about 20 percent over a spatial scale of only 
a few pixels and it does not move relative to the detector. 
Thus, using different profiles for the object, arc and fiat spec- 
tra results in inaccurate flat-fielding and the introduction of 
spurious high-frequency noise into the spectra. This appears 
to be enough to reduce the precision of the radial velocities 
that can be measured from these spectra to 2 - 3kms~^. 

For these reasons we decided to develop our own method 
for extracting the spectra. The key feature of the method 
is to use the flat-field images to create an empirical model 



of the spatial profile for the spectra. This spatial profile can 
then be used to perform a weighted extraction of the spectra 
and to identify pixels affected by cosmic-rays in the images. 
The same weights can be used for the object frames, flat- 
fleld frames and arc frames, so the flat-fielding process does 
not introduce high-frequency noise into the spectra, as is the 
case for optimal extraction with girBLDRS. One disadvan- 
tage of this method is that small shifts in position between 
the flat-field images and the object frames mean that the 
weights applied in the extraction are not quite optimal. This 
results in a small reduction in the signal-to-noise ratio of the 
extracted spectra. A more important effect is that the flux 
at each wavelength is slightly underestimated and that this 
flux-deflcit varies with wavelength and between spectra. For 
our GIRAFFE spectra we find that the flux-deflcit is a few 
percent and that it varies smoothly wavelength. This only 
becomes a problem during the sky-subtraction phase of the 
data reduction. The flrst step in sky subtraction is to form 
an average sky spectrum from those fibres that were pointing 
at blank areas of the sky during the exposure. It is normally 
possible to simply subtract this average sky spectrum from 
the spectra of targets obtained in the same pointing. In our 
case, we first had to calculate an optimum scaling factor for 
each spectrum to be applied to the average sky spectrum 
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prior to subtraction. This was done by finding the scaling 
factor that minimized the root-mean square (RMS) differ- 
ence between the sky-subtracted spectrum and a smoothed 
version of this spectrum. 

We calculated our own dispersion solution for the spec- 
tra from the ThAr spectra obtained on the same day 
as the actual observations. Typical shifts over a 12 hour 
timesc ale for the gir affe spectrograph are less than 0.2 
pixels ijPasauini et al.ll2004ll . We used a 6-th order polyno- 
mial fit to the positions of 22 unsaturated arc lines in the 
arc spectra. The worst RMS residual was 0.092A, the me- 
dian RMS residual was 0.0045A. The mean dispersion of the 
spectra is approximately 4.9 kms~^ pixel^^. 

In order to correct the spectra for telluric absorption 
we used synthetic absorpti on spectra from a 6-layer model 



of the Earth's atmo sphere llNicholk [1988') and the hitran 
molecular database (|Rothman et al. 12005). The parameters 



of the model were optimized by fitting a UVES spectrum ob- 
tained at the same time as the giraffe observations. The 
optimum fit was achieved by minimizing the mean absolute 
deviation from a low-order polynomial fit to the uvES spec- 
trum after dividing through by the model spectrum. The 
fit to the telluric absorption was very good (Fig. [2]). This 
synthetic telluric spectrum was convolved with a Gaussian 
profile to match the resolution of the giraffe spectra. The 
synthetic telluric spectrum was then divided into the target 
spectra. This removed all visible traces of telluric contamina- 
tion. All GIRAFFE spectra were interpolated onto a uniform 
velocity scale of 4.93 kms~^ pixel"^ using 3200 pixels cov- 
ering the wavelength range 8061.1 -8496.6A. We excluded 
spectra with a mean signal-to-noise ratio less than 5 from 
our analysis and also excluded stars with fewer than two 
such spectra. Examples of the resulting spectra for several 
targets are shown in Fig. [3] 



3 ANALYSIS 

Stars ar e identified in this pa per by the J2000 coordinates as 
listed in lKenvon et al.l l|2005h or lBarrado v Navascues et al.l 
l|2003l ) truncated to one decimal place in right ascension and 
truncated to the nearest arcsecond in declination. 



3.1 Radial velocity measurements 

Radial velocities for all targets were measured using cross- 
correlation against a template spectrum of the brown dwarf 
star UScoCTIO 055. This star is a visual binary with two 
similar components separat ed by only 0.12 arc sec and a com- 
bined spectral type of M5.5 (|Kraus et ahlboOSl ). We obtained 
5 UVES spectra of this star from the ESO archive and formed 
the median average spectrum. We re-binned this spectrum 
onto the same wavelength scale as the GIRAFFE spectra and 
convolved this spectrum with a Gaussian with full-width 
at half-maximum (FWHM) of 3 pixels to match approxi- 
mately the resolution of these spectra. Regions of the spec- 
tra affected by strong sky line emission were excluded from 
the cross-correlation. The radial velocity (RV) derived from 
the position of the peak of cross-correlation function (CCF) 
and its error was measured using a parabolic fit to the three 
points at the top of the C CF. The RV of UScoCT IO 055 was 
taken to be —6. 38km s"'^ (jKurosawa et al.ll2006l ). 



Table 2. Radial velocities, Vr, of stars in a Ori and A Ori. The 
date of observation is given as modified heliocentric Julian date 
(MHJD). The standard error of the RV given here includes the 
"external error" and zero-point correction measured from the sky 
lines described in the text. The full-width at half-maximum of 
the cross-correlation function is given in the final column. The 
full version of this table is only available in the on-line version 
of this paper. 



Star 




MHJD 


Vr 




FWHM 








(kms-1) 




(kms-i) 


J053557.0H 


-094652 


53657.309 


28.89 ± 


0.67 


63 


J053557.0H 


-094652 


53691.200 


28.52 ± 


0.60 


63 


J053557.0H 


-094652 


53734.145 


28.74 ± 


0.68 


63 


J053539.4H 


-095032 


53657.309 


27.22 ± 


0.68 


69 


J053539.4H 


-095032 


53691.200 


27.31 ± 


0.61 


70 


J053539.4H 


-095032 


53734.145 


27.82 ± 


0.69 


68 


J053530.4+095034 


53657.309 


28.70 ± 


0.71 


116 


J053530.4H 


-095034 


53691.200 


28.09 ± 


0.66 


116 


J053530.4H 


-095034 


53734.145 


28.77 ± 


0.72 


113 



The precision with which the peak of the CCF can be 
measured is unlikely to be a true refiection of the accuracy 
with which we can measure the radial velocities from our 
spectra. In order to quantify the accuracy of our radial ve- 
locities we measured the radial velocities of the night sky 
emission lines in our spectra. We used the night sky emis- 
sion line spectrum observed with the uvES spectro graph as a 
template in the cross-correlation (|Hanuschikl2003l ). We mea- 
sured all the spectra observed at a given pointing, including 
object spectra prior to sky-subtraction, and calculated the 
mean and standard deviation of the resulting radial veloci- 
ties. The standard deviation, a"sky, is given in Table [T] The 
mean RV shift of the sky-line spectra for a given pointing 
ranges from — 1.28kms~^ to — 0.25kms~^. We have sub- 
tracted this mean radial velocity shift from the measured 
radial velocities of our targets, i.e., the night sky emission 
line spectrum is used to define the zero-point of the stellar 
radial velocities at each pointing. We have also added the 
value of (Jsky for each frame in quadrature to the precision of 
the stellar RV measured from the cross correlation function. 
For 85 percent of our RV measurements o-gky is the dominant 
term in the uncertainty. The resulting RV measurements 
with their standard deviations are given in Tabled 



3.2 Variability criterion for radial velocities 

For each star we have A'rv RV measurements Vr,i each with 
standard error ai. We calculate the weighted mean RV for 
each star, Vr, and then calculate the chi-squared statistic 
for Vr as a model for the observed radial velocities, i.e.. 



JVrv 



1 = 1 



In order to identify stars with variable radial velocities we 
calculate the probability p of observing this value of or 
greater from a sample of normally distributed random ob- 
servations with mean Vr and standard errors ai. Our cri- 
terion for identifying stars with variable radial velocities 
is log(p) < —4. The probability that one or more stars 
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Table 3. Summary of measurements for each target. The full version of this table is only available in the on-line version of this paper. 



Star Ic Afrv Vr log(p) EW(Nal) TiO(8442) 

(mag) (kms"^) (A) 



J053557.0+094652 


14.02 


3 


28.70 


± 


0.37 


-0.04 


2.34 


± 


0.01 


0.677 


± 


0.001 


J053539. 4+095032 


14.06 


3 


27.44 


± 


0.38 


-0.10 


2.34 


± 


0.02 


0.667 


± 


0.002 


J053530.4+095034 


14.10 


7 


28.98 


± 


0.44 


-1.76 


2.33 


± 


0.01 


0.668 


± 


0.001 


J053502. 7+095647 


14.16 


4 


30.06 


± 


1.01 


-4.57 


2.33 


± 


0.02 


0.663 


± 


0.001 


J053408. 4+095125 


14.17 


4 


6.52 


± 


0.34 


-0.05 


2.61 


± 


0.02 


0.637 


± 


0.001 


J053426. 0+095149 


14.36 


4 


-0.65 


± 


0.34 


-0.08 


2.85 


± 


0.01 


0.657 


± 


0.001 


J053555. 6+095053 


14.38 


3 


27.26 


± 


0.39 


-0.08 


2.53 


± 


0.02 


0.682 


± 


0.002 



: J053502. 7 + 095647, log{p)= -4.6" 

: I 


I ^ 




I: 


" J053530. 4 + 095034, log{p)= -1.8 
r I 




L J 


I " 


I i-: 



53680 53700 
MJD 



H 

m 



J- 



Figure 4. Example radial velocity measurements as a function 
of time. The star identifier and value of log(p) are indicated in 
each panel. The weighted mean radial velocity is indicated by a 
dotted line. 



Figure 6. The ratio of the standard errors of the mean radial 
velocities, (Text and cri„t, as a function of Ic magnitude. The width 
of the CCFs is indicated as follows: crosses, FWHM < 80kms~^; 
diamonds, 80kms"^< FWHM < 100 km s"^ asterisks, FWHM 
> 100 km s"'^. Stars with significantly variable radial velocities 
are highlighted (squares). 




: : I : : : I : : : I : : : I : : : 

0.0 0.2 0.4 0.5 0.8 1.0 

P 



Figure 5. Left panel; Cumulative distribution of p. The difference 
between this distribution and a uniform distribution (dashed line) 
is not significant. 



are incorrectly identified as having variable radial veloci- 
ties by chance due to statistical fluctuations (assumed to be 
normally distributed) in our sample of 218 stars is about 
2 percent. The values of A'^rv, log(p) and Vr for each tar- 
get are given in Table [3] There are two ways to calculate 



the standard error of the weighted mean, the external error 
based on the scatter of the data, 

and the internal error based on the standard errors only. 



fint ~ 




The value given in Table |3] is the larger of these two values. 

There are 12 stars in our sample that have variable ra- 
dial velocities according to our criterion log(p) < —4. The 
properties of these stars are listed in Tabled Stars in which 
both components are visible in the spectra or CCFs are iden- 
tified as SB2, those showing only a single spectrum are iden- 
tified as SBl. Examples of radial velocities for stars with 
log(p) < —4 and log(p) > —4 are shown in Fig. |4] 

If our values of ai are good estimates of the true uncer- 
tainty on each RV measurement and the majority of stars in 
our sample are non-variable, then we would expect that the 
distribution of p will be uniform. The cumulative distribu- 
tion of our measured p values is compared to a uniform dis- 
tribution in Fig. [5l We have tested the hypothesis that these 
two distribution are equal using the Kolmogorov-Smirnov 
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Figure 7. The range in full-width at half maximum (FWHM) of 
the cross-correlation function (CCF) for each target as a function 
of Ic magnitude. Stars with variable radial velocities are high- 
lighted (squares). The SB2 binary J054001. 0-021959 discussed 
in the text is marked with an asterisk. 



Figure 9. The distributions of weighted mean radial velocity for 
our targets in A Ori and cr Ori. Dashed lines shows the selection 
criteria for assigning stars to Group 1 or Group 2 for cr Ori and for 
selecting non-members for both clusters, as described in section 



test and find that there is no evidence for any significant 
difference between them. 

We have also tested the reliability of our ai values by 
considering the ratio Z = (Tcxt/o"int. If the values of ai are 
reliable then we expect Z ~ 1. More precisely, for sam- 
ples drawn from normal dis tributions t he moan value of Z is 
1 with standard error l/A/2(iV„ - 1) ( Topping ITqgI) . The 
values of Z for our RV measurements are shown as a function 
of Ic magnitude in Fig. [S] It can be seen that the values of Z 
axe indeed close to 1 (with the exception of the spectroscopic 
binaries, of course) and that there is no significant trend of 
Z with Ic- Stars with broad spectral lines are highlighted 
in this figure so that it can also be seen that the values of 
Z are not significantly different for rapidly rotating stars 
compared to other stars in the sample. 



3.3 Binary stars identified from variable line 
width 

In Fig. [3 we show the range in the FWHM of the CCFs 
against Ic magnitude. We have used this plot to identify 
potential SB2 binary stars in which blending of the spec- 
tra from two similar components results in variations in the 
width of the spectral lines with little change in the radial 
velocity measured from the peak of the CCF. Several of the 
SB2 binaries identified above are recovered by this method. 
The star J054001.0-021959 has a large range in FWHM 
compared to other stars of the same Ic magnitude. The four 
spectra of this star all have high signal-to-noise. There is a 
definite broadening and asymmetry in one spectrum (Fig.jS]). 
In addition, the value of log(p) for this star is close to our 
criterion for variable radial velocities. This star appears to 
be an SB2 binary. One other star (J053522. 5-1-094501) has a 
range m FWHM > 20kms-\ but this is a result of two spec- 
tra of the seven obtained having low signal-to-noise. There 
is no sign of asymmetry in the CCF and the radial veloci- 
ties of this star show no hint of variation (log(p) = —0.09). 
Despite the large range of FWHM measured for this star, 
there is no strong evidence that it an SB2 binary star. 




0-6 0,7 0,8 0,9 

Ti0(a442) 



Figure 10. The value of EW(Nal) versus the TiO(8442) spectral 
index for our targets. The plotting symbols used indicate binarity 
or membership of the a Ori or A Ori clusters or non-membership 
of these clusters based on the radial velocity (if non- variable) , as 
indicated in the legend. Stars below the dashed line are considered 
to be non-members. Error bars are only shown in cases where they 
are larger than the plotting symbol used. 

3.4 Membership criteria 

In this section we describe the criteria we have used to iden- 
tify members of the a Ori and A Ori associations. The prin- 
cipal means of identifying members of these clusters is the 
mean RV of the star, but we have also used the equivalent 
width of the Nal doublet, EW(Nal), and the strength of the 
TiO band at 8442A, TiO(8442), as additional membership 
criteria. 

The distributions of weighted mean radial velocities, 
Vr, are shown separately for stars near cr Ori and A Ori in 
Fig. [9] The bi-modal di stribution of radial v elocities for stars 
near a Ori discussed bv lJeffries et al.l (|2006l ) is apparent. We 
follow the convention in that paper of assigning stars with 
Vr < 27kms~^ to Group 1 and stars with Vr ^ 27kms~^ 
to Group 2. The interpretation of these groups favoured 
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Figure 8. Spectra of the spectroscopic binaries discovered by our survey around the Nal 8200 feature. The spectra have been normahzed 
using a low-order polynomial and offset for clarity. Spectra are labeled by the modified heliocentric JD of observation. The wavelengths 
of the Nal doublet at the mean radial velocity of the star are marked by vertical lines. 
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Figure 8 — continued 



by I Jeffries et alj is that Group 1 are members of either the 
Orion OBla or OBlb association while Group 2 are a sep- 
arate cluster of stars associated with the star a Ori. Group 
2 are concentrated spatially around the star a Ori and have 
similar mean radial velocity to it. Group 2 are younger on 
average than Group 1, although considerable overlap is pos- 
sible. 

For stars near a Ori we identify non-members using the 
criterion Vr > 35km s~^ or Vr < 20km s^^. For stars near 
A Ori we identify non-members using the criterion Vr > 
32km s~^ or Vr < 22km s~^. These limits are indicated in 
Fig- El The application of these criteria to stars with variable 
radial velocities are discussed in more detail in section |4l 

The TiO bands around 8450A are good indicators of 
effective temperature for M-dwarfs in the sense that they in- 



crease in strength for cooler stars (|Mohantv et al.ll2004h and 
are insensitive to reddening. We have measured the strength 
of the band at 8442A in our spectra using the ratio of the 
counts detected in the wavelength ranges 8437-8442A and 
8442 - 8450A. The measurements were made on the median 
average spectrum of each star. We denote the value of this 
ratio as TiO (8442). 

The equivalent width of the Nal doublet at 8190A 
is sensitive to the surface gravity in M-type stars 
jSchiavon et al.iri997l ) in that stars with higher surface grav- 
ities show stronger NaT absorption. The surface gravity of 
members of the a Ori and A Ori clusters is expected to 
be log g = 3-4, whereas a typical M-type giant will have 
log g ~ 1. Thus, the equivalent width of the Nal doublet, 
EW(Nal), can be used to identify background giants in our 
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Table 4. Properties of the binaries discovered by our survey. The standard error in the mean for the measured radial velocities 
is listed under arv Non-members of the clusters are indicated in the final column by the note 'nm'. 



Star Ic Rc-Ic Nrv Vr (Trv log(p) EW(Na I) TiO(8442) Notes 

(mag) (mag) (kms""*^) (kms"'^) (A) 



J053502. 7+095649 


14.16 


1.30 


4 


30.06 


1.01 




-4.6 


2.33 ± 0.02 


0.663 ± 0.001 


SB2 


J053456.3+095503 


14.54 


1.36 


4 


26.20 


2.32 




28.8 


2.21 ± 0.01 


0.656 ± 0.001 


SBl 


J053612. 1+100056 


14.60 


1.31 


3 


28.66 


2.68 




-8.9 


2.23 ± 0.10 


0.666 ± 0.008 


SBl 


J053443. 9+094835 


15.20 


1.69 


4 


22.56 


3.48 


< 


-42 


2.66 ± 0.02 


0.732 ± 0.002 


SB2 


J053455. 2+100034 


15.23 


1.72 


4 


27.92 


2.20 




15.6 


2.37 ± 0.02 


0.710 ± 0.002 


SB1_^ 


J053845. 6-021157 


14.48 


1.18 


4 


4.11 


3.52 


< 


-42 


1.99 ± 0.02 


0.633 ± 0.001 


SB1~ 


J053801. 0-024537 


14.46 


1.57 


4 


46.12 


9.20 


< 


-42 


2.09 ± 0.01 


0.653 ± 0.001 


SB2 


J054001. 0-021959 


15.02 


1.90 


4 


30.80 


0.90 




-3.5 


2.38 ± 0.02 


0.726 ± 0.003 


SB2_^ 


J053823. 5-024131 


15.15 


1.69 


4 


30.23 


4.26 


< 


-42 


2.71 ± 0.02 


0.690 ± 0.002 


SB2~ 


J054052. 5-021653 


14.30 


1.13 


4 


152.28 


8.25 


< 


-42 


0.52 ± 0.02 


0.628 ± 0.001 


SBl, nm 


J053743. 5-020905 


14.48 


1.36 


4 


51.02 


4.76 


< 


-42 


2.65 ± 0.02 


0.630 ± 0.001 


SB2 


J053838. 1-023202 


16.59 


1.73 


4 


30.79 


2.65 




19.5 


2.05 ± 0.07 


0.659 ± 0.006 


SBl 


J053805. 6-024019 


14.13 


1.21 


4 


7.68 


2.53 




41.1 


2.43 ± 0.01 


0.644 ± 0.001 


SBl 



" Hint of companion in CCF. 

Identified from variation in width of CCF. 
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Figure 11. The TiO(8442) spectral index versus I-band magni- 
tude for our targets. The plotting symbols used indicate binarity 
or membership of the a Ori or A Ori clusters or non-membership 
of these clusters based on the radial velocity or spectral type (see 
Fig, lion , as indicated in the legend. Error bars are only shown in 
cases where they are larger than the plotting symbol used. 

sample. It is also possible to identify contamination of the 
sample by dwarf stars with log g ^ 4.5. We measured the 
value of EW(Na I) from the median average spectrum of each 
star by numerically integrating the area under our spectra 
in two regions ± 120km s~^ wide around the centre of each 
Nal line after normalizing the spectra by the clipped mean 
value of the spectrum in the region 8188 -8192 A. 

The value of TiO(8442) is plotted against EW(Nal) in 
Fig. 1101 There is a clear division between the bulk of our 
targets and the small group of stars with low EW(Nal) 
and TiO(8442) values. Most of these stars are non-members 
based on their mean RV. We therefore identify all stars 
below the dashed line in this figure as "spectral-type non- 
members" . 

From Figs. 1111 and [T2l we see that very few stars that 



Figure 12. The equivalent width of the Nal doublet for our 
targets as a function of the I-band magnitude. The plotting sym- 
bols used indicate binarity or membership of the a Ori or A Ori 
clusters or non-membership of these clusters based on the radial 
velocity or spectral type (see Fie. HOC , as indicated in the legend. 
Error bars are only shown in cases where they are larger than the 
plotting symbol used. 

satisfy the RV and spectral type criteria for membership 
have discrepant values of EW(Nal) or TiO(8442) for their 
Ic magnitude. There are two faint stars in the sample that 
have large values of TiO(8442), but there is a large scatter 
in this index among the faint stars in our sample so we do 
not consider this to be a reason to exclude these stars as 
members of the a Ori cluster. 

The binary star with a discrepant value of EW(Nal) is 
J054052. 5-021653. This SBl binary star has weak, narrow 
Nal absorption lines characteristic of a giant star. The mean 
RV has a rather large uncertainty but is also clearly incon- 
sistent with membership of the a Ori cluster. We conclude 
that this is a background giant star. 

The spectra of the star J054034.5-020606 extracted us- 
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Figure 13. The I-band mass-magnit ude relations interpolated 
from tlie models of Bara ffe et alj l ll998l) for the following distances 
and ages: d=330pc, agc=10Myr [a Ori 1); d=440pc, age=3Myr 
(o- Ori 2); d=400pc, age=5 Myr (A Ori). 

ing our empirical weighting scheme do not satisfy our crite- 
rion for inclusion in the sample because the signal-to-noise 
is less than 5 but we did notice that this is an SB2 spec- 
troscopic binary from an analysis of the spectra extracted 
using girBLDRS. These spectra are shown in Fig. [S] This 
is the faintest binary detected in our survey (Ic= 18.38) so 
the spectra are quite noisy, but the SB2 nature of this star 
can be seen and is very obvious in the CCF. The values 
EW(NaI)= 3.42 ± 0.23 of TiO(8442)= 0.593 ± 0.026 shows 
that this star has the wrong spectral type to be a member 
of the a Ori cluster. We have also inspected the position of 
this star in the V v. V— Ic col our-magnitud e diagr am of the 
0" Ori region using the data of lMavne et al.l l|2007l ), where it 
is clearly below the sequence of cluster members. We con- 
clude that this is a background dwarf binary star for which 
the combination of spectral type, distance and reddening 
place the star within the band of cluster members in the 
Ic V. Rc — Ic colour-magnitude diagram. 

Having applied these membership selection criteria we 
find that in our sample there are 64 members of A Ori, 34 
members of a Ori Group 1 and 86 members of a Ori Group 
2, excluding stars with variable radial velocities that are 
discussed separately below. 



4 BINARY STARS IN OUR SAMPLE 

In this section we discuss the stars in our sample with 
variable radial velocities or variable line profiles that are 
members of the a Ori or A Ori clusters. The values of Vr, 
TiO(8442), etc. for these stars are given in Table ID The 
spectra of these stars in the region of the Nal doublet are 
shown in Fig. |S] 

We also note in Table |4] whether only one set of spectral 
lines are visible (SBl) or that the spectra indicate the pres- 
ence of two stars (SB2). A star may be an SBl binary either 
because the companion is too faint to be detectable or be- 
cause the lines of the stars are blended, or both. We created 
a set of synthetic binary star spectra using combinations of 
observed single star spectra over a range of luminosity ra- 
tio and radial velocity difference to estimate the lower limit 



to the magnitude difference in the SBl stars. If we assume 
that neither star is rapidly rotating then we estimate that 
any companion to these stars is more than 2.5 magnitudes 
fainter than the primary. From the magnitude-mass relation 
shown in Fig. [13] we estimate that the corresponding limit 
to the mass ratio is q ^ 0.25. 

In general, the IcRc photometry, the values of 
TiO(8442), EW(Nal) and the mean RV of these stars are 
all consistent with membership of either the A Ori cluster 
or one of the a Ori clusters. We discuss exceptions to this 
general rule and any other points of interest for these stars 
below. 

J053502.7-H095649 iBarrado v Navascues et"al] (|2004| ) 

consider this star (LOri-CFHT 043) to be a member of 
the A Ori cluster on the basis of the available RcIcJHK 
photometry. 

J053443. 9+094835 iBarrado v Navascues et al.1 (|2004l ) 
note that the position of this star (LOri-CFHT 069) in 
the Iq v. Ic— Ks colour-magnitude diagram is inconsistent 
with cluster membership. However, this is simply a result 
of the contribution of two stars of similar brightness to the 
flux at I-band in this SB2 binary (Fig. [Sjl. Correcting the 
I-band magnitude by 0.75 magnitudes places this star in a 
position entirely consistent with cluster membership in this 
colour-magnitude diagram. 

J053455.2-f 100035 .Barrado v Navascues et al.1 (|2004| ) 
note the presence of Ha emission in their low resolution 
spectrum of this star (LOri-CFHT 075). The position of 
this star in the Iq— Ks v. H— Ks colour-colour diagram 
and the Ic v. Ic— Ks colour-magnitude diagram are also 
consistent with cluster membership and with the spectra l 
type of M5 assigned by IBarrado v Navascues et al.l (|2004l ) 
so it is unclear why they note it as being a non-member on 
the basis of this information in their Table 2. The spectral 
lines for this star show rotational broadening. We compared 
the spectra of this star to those of a narrow-lined star of 
similar spectral type to which we had applied a rotational 
broadening function for various values of the projected rota- 
tional velocity, Vrot sini. From this comparison we estimate 
that the projected rotational velocity of J053455. 2-1-100035 
is Vrot sin i ~ 65km s~^. Only one set of spectral lines are 
visible in our spectra but there is an asymmetry in the CCF 
in the form of a blue-wing, particularly when the measured 
RV corresponds to a red-shift. The suggests that the fainter 
component in this binary is detected but unresolved in the 
I-band spectra. 

J053845. 6— 021157 The range in radial velocities we 
have observed for this star is —6.9 to 11.2km s^^. This is 
consistent with this star being a member of a Ori Group 1 
if it has a semi- amplitude « 30km s~^. Members of Group 

1 are approximately twice as common as members of Group 

2 in the 25arcmin giraffe field used for the observations 
of this star. The orbital period is required to be rather 
short (P ^ 1 day) to reconcile the observed radial veloci- 
ties of this star with a mean velocity consistent with cluster 
membership and the expected mass ratio for an SBl binary 
(g^O.25). 

J053801. 0-024537 The mean RV estimated from the 
spectra is consistent with membership of either a Ori Group 
1 or cr Ori Group 2. There are approximately 10 times as 
many members of Group 2 as Group 1 in this field close to 



12 P.F.L. Maxted et al. 



a Ori so we conclude that it is approximately 10 times more 
likely that this star is a member of Group 2 than Group 1. 
The values of EW(Nal) and TiO(8442) for this star were 
measured from the spectrum observed near conjunction. 

J054001. 0-021959 iKenvon et al.l (|2005h detected strong 
lithium absorption in the spectrum of this star (KJN2005 6), 
which indicates that this star is younger than 20Myr, as 
expected for a member of the a Ori cluster. The mean RV 
of this star suggests it is a more likely to be a member of 
(J Ori Group 2 than Group 1. 

J053823. 5-024131 The mean RV measured from the 
three spectra in which the two components are unresolved 
is 26.5 ± 0.9km s^^, which is close to the dividing line at 
27kms~^ between Group 1 and Group 2. The proximity of 
this star to a Ori makes it approximately 10 times more 
likely that this is a member of Group 2 than Group 1. This 
star (BMZ2001 S Ori J053823.6-024132) was listed as a can- 
didate member of the a O r i clus ter on the basis of IJHK 
photometry bv'B eiar et al.l (|2004l '). 

J053743. 5-020905 We measured the value of EW(Nal) 
and TiO(8442) for this SB2 binary from the spectrum taken 
near conjunction, so the values represent an average value 
for the two stars, which are of comparable brightness. The 
RV measured from this spectrum is 42kms~^, which is out- 
side the range we have defined for membership of the a Ori 
cluster, but it is difficult to establish whether this is an ac- 
curate estimate of the mean RV of this star. Further ob- 
servations will be required to establish the true mean RV 
of this binary star in order to check that it is consistent 
with this star being a cluster member. Members of Group 1 
are approximately tw ice as common in this field as Group 
2 (| Jeffries et al.|[2006l ). For the purposes of this paper we 
assume that this star is a cluster member. 

J053838. 1-023202 The mean RV of this SBl binary and 
its proximity to a Ori suggest that it is a member of Group 
2. 

J053805. 6-024019 The range of radial velocities ob- 
served in this SBl binary star is 4.4 to 14.5km s"'^, which is 
consistent with membership of a Ori if the semi-amplitude 
of the spectroscopic orbit is if ~ 20km s~^. The proximity 
of this star to a Ori makes it approximately 10 times more 
likely that this is a member of Group 2 than Group 1. 



In summary, we have identified 5 spectroscopic binary 
members of the A Ori cluster and 7 spectroscopic binary 
members of the a Ori cluster. Among the a Ori spectro- 
scopic binaries, 5 stars are likely to be members of Group 
2, 1 is likely to be a member of Group 1 and 1 star cannot 
be assigned to either group. Of the 12 spectroscopic binary 
cluster members, 6 are SB2 binaries with stars of compa- 
rable brightness in the I-band. There is 1 SBl binary with 
broad lines due to rotation so it is not clear whether this is a 
genuine SBl binary or an unresolved SB2 binary, although 
there is a hint of the companion in the CCF for this star. 
There are 3 SBl binaries with narrow lines for which we can 
say the companion is likely to be more than 2.5 magnitudes 
fainter than the primary so that the mass ratio is less than 
about 0.25. We have also identified two spectroscopic bina- 
ries that are not members of either a Ori or A Ori. 



4.1 The distribution of binaries with magnitude 

It is notable that all 12 spectroscopic binary stars we have 
identified among the 196 cluster members are brighter than 
Ic ~ 16.6, and 11 are brighter than Ic ~ 15.25. There are 
68 cluster members brighter than Ic = 15.25 including these 
binary stars. At face value, this suggests a close binary frac- 
tion of about 16 percent above this limit and less than a few 
percent below this limit. Of course, the sensitivity of our 
survey decreases for fainter stars because the signal-to-noise 
of the spectra is less and the orbital speeds decrease with 
mass for a given semi-major axis or orbital period. 

We have used the RV data for the 11 binary clus- 
ter members we detected from their RV variations to es- 
timate the sensitivity of our survey to binaries of this type 
as a function of I-band magnitude, Pompiricai. We exclude 
J054001. 0—021959 from this analysis because it is unclear 
how to simulate the method by which we detected this bi- 
nary, i.e., from the variation of the width of its CCF. For 
every combination of variable and non-variable star in our 
survey we have created a synthetic RV dataset using the RV 
errors in the non-variable star dataset and the radial veloci- 
ties of the binary star scaled as follows. We first subtract our 
best estimate of the systemic radial velocity for the binary. 
We then estimate the total mass of the binary, tut assuming 
that the stars are identical for the SB2 binaries or that the 
companion is 2.5 magnitudes fainter than the primary for 
the S Bl binaries. The mas ses are estimated using the mod- 
els of iBaraffe et al.l (|l998l ) based on the Ic magnitude of 
the target. The relations between Ic and mass are shown in 
Fig. [T31 We have assumed the following distances and ages 
for the clusters: d=330pc, age=10Myr {a Ori 1); d=440pc, 
age=3Myr (a Ori 2); d=400pc, age=5Myr (A Ori). We 
then repeat the calculation to find the total mass, Mt, of 
a similar binary star with the same Ic magnitude as the 
si ngle star , ms- We then multiply the radial velocities by 
•y/ ms/rriT- This is equivalent to assuming that the distri- 
bution of semi-major axis is independent of mass. We then 
apply the same detection criterion as before (log(p) < —4) 
to the synthetic datasets. In cases where A^'rv is less for the 
single star than the binary star we use the average number 
of detections for all combinations of A^rv synthetic radial 
velocities from the combinations available. For cases where 
Arv is larger for the single star than the binary star we use 
the Arv measurements with the lowest RV errors to create 
the synthetic dataset. The number of synthetic RV data sets 
which satisfy our variability criterion then gives an estimate 
of the detection efficiency for each star, e.g., if 6 of the 11 
synthetic RV data sets for a star satisfy our variability cri- 
terion, then the same observations of a binary star similar 
to those discovered in this survey with the same Ic magni- 
tude would, on average, have detected a significant RV shift 
about 55 percent of the time. The values of Pcmpiricai calcu- 
lated in this way is shown in Fig. [14] as a function of the Ic 
magnitude. 

The normalized cumulative distribution function for 
these detection efficiencies is shown as a function of Ic mag- 
nitude in Fig. 1151 Also shown is the cumulative distribu- 
tion function for the Ic magnitude of the binary cluster 
members excluding J054001. 0—021959. If the binary frac- 
tion and semi-major axis distribution for binaries is inde- 
pendent of mass then these two distribution should be the 
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same. It is clear that the two distributions are not the same. 
The Kolmogorov-Smirnov test apphed to these distributions 
gives a 99.7 percent significance level to the difference in 
these distributions. There is a chance that a few of the stars 
we have identified as SBl binary stars are the result of spu- 
rious RV shifts due to errors in the analysis or instrumental 
effects or the result of intrinsic variability of a single star. 
If we take a very cautious approach and repeat this anal- 
ysis using the detections of SB2 binaries only we find that 
Kolmogorov-Smirnov test gives significance level of 97.6 per- 
cent. 



4.2 The binary fraction 

We used a Monte-Carlo simulation to calculate probability 
distribution for the binary fraction given the 11 binaries we 
have discovered from their variable radial velocities and the 
detection efficiency for each star, Pompiricai, calculated above. 

The results are shown in Fig.[T6]for the whole sample, a 
'bright' sample (Ic < 16.9) and a 'faint' sample (Ic ^ 16.9). 
The division here between the bright and faint samples 
has been chosen to correspond to the widely accepted di- 
vision between low mass star s and VLM stars at a mass of 
O.lM0 l|Burgasser et al.|[2007l '). The mean value of Pempiricai 
for the 145 stars in the bright sample is 0.89. For the 51 
stars in the faint sample the mean value of Pcmpiricai is 0.53. 
The spectroscopic binary fraction for the bright sample is 
9.5^^2 g percent. The 90 percent confidence upper limit to 
the spectroscopic binary fraction for the faint sample, given 
the assumptions above, is 7.5 percent. These figures apply to 
binary stars of the type discovered by our survey only, i.e,. 
we have not attempted to apply a correction to this spec- 
troscopic binary fraction for the binaries at longer orbital 
periods that are not detected by our survey. The hypoth- 
esis that /bright and /faint are equal can be rejected with 
90 percent confidence. The value of Ic ~ 16.9 that is used 
as the dividing line between the faint and bright samples 
is arbitrary. It would be possible to achieve a higher level 
of significance by setting the limit between our bright and 
faint samples closer to magnitude of the faintest binary in 
our sample (Ic = 16.5), but this would be an equally arbi- 
trary value and so it would be hard to justify the apparent 
increase in statistical significance. 

4.3 The detection efficiency as a function of 
binary separation 

We have used another Monte Carlo simulation to estimate 
the range of semi-major axis, a, over which binary stars can 
be detected by our survey. 

If binarity is the only cause of variable RVs, the prob- 
ability that a given target is flagged as an RV variable is 
given by £(,Pdctcct + (1 — eb)10~*, where ei, is the overall bi- 
nary fraction and pdotcct is the probability that logp < —4 
for the object assuming that it is a binary. 

We have used a Monte Carlo simulation to calculate 
the value of Pdctcct as a function of semi-major axis, a, for 
every star in our sample given various assumptions about 
the distribution of binary properties. The simulation gener- 
ates 65536 virtual binaries and predicts the RV of the more 
massive component at the same times of observation as the 
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Figure 14. The detection efficiency as a function of magnitude 
for our survey based on the radial velocities of the binaries iden- 
tified from their radial velocity variations. 




Figure 15. The normalized cumulative distribution of the detec- 
tion efficiency of our sample as a function of Ic magnitude based 
on the measured radial velocities of all binaries detected (solid 
line) or the SB2 binaries only (dashed-dotted line) compared to 
the normalized cumulative distribution of Ic magnitude for all 
the binaries discovered by our survey (dashed line) and the SB2 
binaries discovered by our survey (dotted line). 

actual observations. The eccentricity, e, mass ratio, q, and 
other properties of the binary star are randomly selected 
from the following distributions. 

Mass ratio, q We have used a 'flat' distribution which 
is uniform in the range q = 0.2-1. We did not con- 
sider the peaked mass ra tio distribution we investigated in 
iMaxted fc JeflFriesI (|2005[ l to be appropriate for the binaries 
we have discovered in this survey because that distribution 
is zero for q < 0.7 whereas some of the SBl binaries we 
have found must have mass ratios q ^ 0.25. The value of 
the mass ratio makes little difference to the value of pdctoct 
in this range. 

Eccentricity, e We have assumed that all bina- 
ries with periods less t han 10 d have circular orbits 
(|Meibom fc Mathieull2005l ). Above this period, we assume 
that the value of e is uniformly distributed in the range 
e = 0-emax where emax = 0.6. We have also performed one 
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Figure 16. The probability distribution for the binary fraction 
in our whole sample (dashed line), the 'bright' sample (solid lines, 
mass > O.IMq) and the 'faint' sample (dashed-dotted line, mass 
< O.IMq). 



set of simulations with emax = in order to investigate the 
effect of assuming that all orbits are circular. 

Primary mass, m We have calculated the mass of the 
primary star based on its Ic magnitude and the models of 
iBaraffe et al] (|l998l ). The mass-magnitude relations shown 
in Fig. [T3] are used to find a primary star mass consis- 
tent with the observed magnitude and the mass ratio of 
the synthetic binary. We allow for an assumed error of 
0.03 magnitudes in the observed Ic magnitude. We find that 
the choice of mass-magnitude relation for each star has a 
negligible affect on our results. For simplicity, we present 
the results assuming a distance of 330 parsec and an age of 
10 Myr for all stars. 

Orbital phase The orbital phase of the binary at the date 
of the first observation is randomly selected from a uniform 
distribution in the range to 1. 

Longitude of periastron, u: For eccentric binaries, uj is 
selected from a uniform distribution in the range to 27r. 

We have calibrated the extent to which blending be- 
tween the components reduces the apparent amplitude of 
the RV variation in a spectroscopic binary. We selected 6 
stars with a range of Ic magnitudes that had typical spec- 
tra for stars of that magnitude. For each pair of stars we 
created simulated binary star spectra in which the single 
star spectra were combined with the appropriate flux ratio 
for their magnitude difference and a range of velocity offsets 
between the stars. We then measured the radial velocity of 
the brighter star in the spectrum in the same way as we did 
for our actual observations. We used these results to cali- 
brate the difference between the true RV and the apparent 
RV of the brighter star caused by blending with the spec- 
trum of the fainter star. We used interpolation within the 
resulting table to adjust the RV of the more massive star in 
each simulated binary star to account for this blending. 

The radial velocities predicted by each trial of the sim- 
ulation are each perturbed by a random value from a Gaus- 
sian distribution with the same standard deviation as the 
random error of the actual observations. We then estimated 
the range of inclinations over which the binary would be 
detected using the same criterion that we applied to our 



Figure 17. The average detection efficiency for our survey as 
a function of semi-major axis (a) for the bright sample (upper 
curves) and the faint sample (lower curves). The solid line cor- 
responds to Cniax = 0.6, the dashed line to emax = and the 
dashed-dotted line shows the effect of neglecting blending in the 
case emax = 0.6 for the bright sample. 



actual data. In the absence of blending this is a trivial cal- 
culation. In the presence of blending the value of may not 
be a monotonic function of inclination, i. We approximate 
the true shape of the relation between and sini using a 
parabolic flt to the values at sini = 0, 0.5 and 1 and use this 
parabolic fit to estimate the range of sin i values over which 
the binary would be detected. 

The average value of Pdctcct calculated in this way for 
all the stars in the bright and faint samples are shown as 
a function of log a in Fig. [TT] From that figure we see that 
assuming circular orbits does not make a large difference to 
the value of Pdetect, but that neglecting blending can lead 
to an overestimate of the detection efficiency by as much as 
15 percent. We can also see that the detection efficiencies 
calculated above using the RV data of the binaries them- 
selves gives a result consistent with the values of pdotcct in 
Fig. [T7] if the semi-major axes of these binaries are in the 
range —2.5 ^ log(a/a-!i) ^ —1.0. The corresponding range in 
orbital periods is 0.1 d ^ P ^ 30 d, which is in good agree- 
ment with the likely orbital periods of the binaries we have 
detected. 



4.4 Comparison with iKenvon et al.l (|2005l ') 

We combined our radial velocity data with those of 
iKenvon et all l|2005l ) to see if the combined data sets would 
yield any further spectroscopic binaries. We did not find 
any new binaries among the 45 stars in common between 
the two surveys. This is, perhaps, not surprising given the 
m uch higher radi al velocity accuracy of our data compared 
to IKenvon et al.l . This is as a result of the higher signal- 
to-noise, higher resolution and superior telluric subtraction 
achievable with giraffe spectrograph compared with the 
WYFFOS spectra available to Kenyon et al. The mean differ- 
ence in the measured radial velocity between the two sets of 
data is 0.72 ± 0.44kms"\ 

The star KJN2005 72 (J053739. 6-0218 26) was iden- 
tified as a possible spectroscopic binary by IKenvon et al.l 
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l|2005l ) based on a RV shift of 35 ± 4kms ^ between two 
spectra obtained on consecutive nights. The radial velocities 
for this star measured from our data are constant to within 
a 2 kms~^ over a baseline of 63 d ays. We conclu d e that the 
radial velocity shift measured by iKenvon et all (|2005l ) for 
this star is lik ely to be spuri ous. 

Similarly. iKenvon et al.l claim a RV shift of 13±6 km s~^ 
for the star KJN2005 74 (J053926. 8-023656) between two 
spectra obtained on consecutive nights. We find the RV for 
this star is constant to within 2km s~^ from 3 spectra with 
a baseli ne of 61 da ys. 

Ke nyon et al] note that the star KJN2005 46 
(J054000. 1—025159) appears to be a member based on the 
presence of the Li I 6707A line in the spectrum and the 
equivalent width of the Nal doublet, but the RV they mea- 
sure (17±3kms~^) is inconsistent with cluster membership. 
They suggest that this may be due to this star being spec- 
troscopic binary. However, we find that this star has a mean 
radial velocity of 30.5km s~^ which is consistent with clus- 
ter membership and is constant to within 1km s~^ from 3 
spectra with a baseline of 61 days. 

In su mmary, it app ears that the radial velocities mea- 
sured by iKenyon et al.l sometimes show spurious shifts of 
approximately 10km s~^. 



5 DISCUSSION 

There are several examples of V LMS and BDs th at are 
clearly spectroscopic binaries. iStassun et al.l (|2006l ) mea- 
sured accurate masses and radii for the brown dwarf pair 
2MASS J05352184-0546085 which is an eclipsing spectro- 
scopic binary with an orbital period of 9.8 days and a total 
mass of O.O88M0 in the Orion Nebula cluster (ONC) . PPl 15 
in the Pleiades was the first brown dwarf confirmed by the 
detection of lithium and is an SB2 spectroscopic binary con- 
taining two brown dwarfs with masses of 60-70 Jupiter 
masses , an orbital period o f 5.8 days and an eccentricity 
of 0.4 (|Basri fc Martini Il999l ). There is no large RV survey 
for spectroscopic binary VLMS and BDs in the Pleiades so 
we do not know whether PPl 15 is representative of the 
binary population in this cluster. Similar arguments apply 
to 2MASS J05352184-0546085 and the ONC. Intriguingly, 
there is a very well populated binary se quence in the colour - 
magnitude diagram for the Pleiades (|Lodieu et al.l [20071 ). 
which suggests a binary frequency of 28-44 percent in the 
0.075 -O.O3OM0 mass range. 

Surveys for spectroscopic binaries among late-M and T 
dwarfs in the field and in nearby clusters have discovered 
a few other SB2 binaries an d several stars that sho w RV 
shifts of a few kms^^ or less. iBasri fc ReinersI l|2006l ) sum- 
marize the results of these surveys and present the results 
of their own survey of 53 VLMS and BDs. From their own 
sample they estimate a binary frequency of 11 percent in 
the separation range 0-6 an. This binary frequency may be 
consistent with our results for VLM binaries (binary frac- 
tion < 7 percent for mass ^ O.IMq) if the distribution of 
a is not strongly biased towards small a. T his is reasonable 
given that none of the binaries detected by lBasri fc Reiner3 
were SB2 binaries and that the typical radial velocity shifts 
detected were small (few kms~^). 

It is much harder to interpret the results of previous 



surveys summarized bv lBasri fc ReinersI (|2006l ). not only be- 
cause of the small number of stars observed but also because 
there is a bias in these surveys due to preferentially observ- 
ing brighter stars or imposing a magnitude limit. This tends 
to increase the number of binaries in these samples, particu- 
larly SB2 binaries. Our survey is much less strongly affected 
by this type of bias. Some SB2 stars will be missing from 
our sample because they exceed our bright magnitude cut- 
off. The lack of faint binaries in our sample means there is 
no compensating gain in SB2 binaries at the faint end of the 
sample. This bias can only increase the statistical signifi- 
cance of the change in binary properties near Ic ~ 16.9 that 
we have discovered. It is difficult to estimate the compen- 
sating bias introduced by excluding stars below the cluster 
sequence in the Ic v. Rc — Ic colour magnitude diagram but 
this effect is expected to be small. 

It is also harder to characterize the sources of noise 
in small surveys. T his is apparent from the results of 
IKenvon et al.l (|2005l ) in which the few RV shifts they mea- 
sured of about 10km s~^ wer e ascribed to binary mot ion, 
both by IKenvon et all and bv lMaxted fc Jeffrie j (120051 '). In 
that case it appears that there is some problem with obtain- 
ing reliable RV measurements from spectra affected by tel- 
luric absorption with low resolution and low signal-to-noise. 
A similar problem applies to the interpretation of radial ve- 
locity shifts ^ 1 kms~^, even for spectra of the highest qual- 
ity. At this level, line profile variations due to chromospheric 
activity and star spots ("jitter") may become important, but 
the extent of this effect has not been well characterized. Our 
survey may also suffer from this problem in the case of some 
of the SBl binaries we have detected. For stars with narrow 
lines and large RV shifts such as J053 805.6-024019 this is 
unlikely to be an issue. For example, Ijoergens fc Mueller! 
(|2003) have shown that the star Cha Ha 8 is a binary star 
with an eccentric orbit (e = 0.49) and a low mass com- 
panion in a wide orbit (a ~ 1 an). The semi-amplitude 
of the orbit is low {K = 1.6km s~^) but can be measured 
reliably because the level of jitter in this young, narrow- 
lined VLM star is much less than lkms~^. It is less cer- 
tain that stars like J053455. 2-1-100035 with broad spectral 
lines and small RV shifts are genuine spectroscopic bina- 
ries. One mitigating factor in our survey is that there are 
several other stars in our survey with broad lines that are 
not detected as binary stars, so it does appear that we can 
measure radial velocities accurate to about 1km s~^ in this 
type of star. Nevertheless, it remains to be established that 
any of the single-lined stars showing RV shifts from this sur- 
vey, the survey bvlBasri fc ReinerSi (|2006) or the survey by 
iGuenther fc Wuchterll 1200^ " are genuine SBl spectroscopic 
binaries with Keplerian orbits. 

Note, however, that the statistical significance of the 
change in binary properties for very low mass stars in this 
survey we have noted remains high (> 97 percent) even if 
we consider only the SB2 binaries. 

It is harder to detect binary stars among rapidly rotat- 
ing stars because the RV measurements from broader CCFs 
is less precise than for narrow CCFs. This does not affect 
our conclusion that there is a change in binary properties 
for very low mass stars compared to low-mass stars because 
rapidly rotating stars are evenly distributed in magnitude 
within our sample and the standard errors of the radial ve- 
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locities we assign to these stars accurately accounts for the 
effects of hne broadening. 

A rapid change in close binary fraction at a given mag- 
nitude is only approximation to what is, in reality, likely to 
be a gradual change in binary properties with mass. We do 
not have enough data to be able to say with any accuracy 
at which mass this change occurs or whether the change oc- 
curs over a small or large range of masses. This is why we 
have presented results with the sample divided into bright 
and faint sub-samples at Ic = 16.9, a magnitude that cor- 
responds to the commonly accepted dividing line between 
VLMS/BDs and low-mass stars at Q.IMq. 

With only 11 binaries and 3 or 4 RV measurements 
per star we cannot say a great deal about the distribution 
of the binaries properties in low mass star with the data 
available so far. However, we do have estimates of the mass 
ratio in these binaries. The mass ratio is clearly g « 1 in 
the case of the SB2 binaries with nearly equal components 
("twins"). For the stars where no companion is visible and 
the spectral lines are narrow we estimate that the companion 
is likely to be at least 2.5 magnitudes fainter than the pri- 
mary and so g ^ 0.25. The star J053455.2-H00035 lies some- 
where between these extremes. The overall picture, then, is 
of a peak in the mass ratio distribution near q = 1 due to 
twins over-laid on a broader distribution increasing towards 
small values of q. This is qualitatively similar to the mass 
ratio distribution foun d for low mass Population I stars by 
iGoldberg etlll l|2003l ) and the mass ratio distribution for 
solar-type stars measured bv lHalbwachs et al.l (|2003l ) for so- 
lar t ype stars in the so lar neighbourhood, Pleiades and Prae- 
sepe. iGoldberg et al.l find that the peak in the q distribution 
due to twins is not present for P op I binaries mor e massive 
than 0.67i\'fQ or for halo stars. iHalbwachs et al] find that 
the peak for q > 0.8 gradually decreases when long-period 
binaries are considered. 

We have estimated the number of binaries we would 
detect if the stars brighter than Ic=16.9 (mass « O.IMq) 
in our sample have the same binary frequency and or- 
bital period distribution as nearby M-dwarfs. We used the 
detection efficiency for each star as a function of log P 
convolved with the Gaussi an distribution for log P from 
iDuguennov fc Mavoil (|l99ll ) toget her with the binary frac - 
tion for M-dwarfs of 42 percent from lFischer fc Marcvl (|l992l ) 
to find that we would have detected 6.5 spectroscopic bi- 
naries an average in our survey given these assumptions. 
At face value, this suggests that the frequency of short pe- 
riod binaries among low mass stars in A Ori and a Ori is 
slightly larger to that for nearby M-dwarfs. However, the 
assumption that the distribution of log P established by 
IDuguennov fc Mavoil for solar-type s tars is appropriat e for 
low mass stars may not be correct. iFischer fc Marcvl find 
that the period distribution for M-dwarfs binaries is simi- 
lar to that for solar type stars but there is only one short 
period M-dwarf binary in the sample so the period dis- 
tribution is poorly ch aracterize d at th e short period end. 
IDuguennov fc Mavoil note that iGriffinI (|l985l ) measured a 
much higher frequency of short period binaries for solar 
type stars in the Hyades than is seen in their field star 
sample. The sam e effect was seen to a lesser extent by 
IHalbwachs et all so both environment and primary mass 
may be important factors in determining the frequency of 
short period binaries in our sample. We do not have suffi- 



cient data in our survey to measure the orbital periods and 
semi- amplitudes of the spectroscopic binaries we have de- 
tected. 

Given the uncertainty in the orbital period distribution 
of low mass stars it would clearly be worthwhile to obtain 
complete orbits for the binary stars we have identified. The 
maximum separations of the lines observed in the SB2 bina- 
ries and the timescale and amplitude of the RV variations 
for the SBl binary stars suggest that the orbital periods 
of these binaries are likely to be in the range from several 
hours to several days. An anonymous referee has raised the 
possibility that the difference in binary fraction we have ob- 
served in the bright and faint samples is due to a distribu- 
tion of log a strongly biased towards \og{a/au) ^ —0.5. It 
can be seen from Fig. [T7] that we have almost no sensitiv- 
ity to binaries within the faint sample but good sensitivity 
to binaries within the bright sample in this log a range. In 
principle, it may be possible to find a distribution of log a 
that can explain the numbers of binaries detected in the 
bright and faint samples for a single value of the binary 
fraction. However, such a distribution would be incompati- 
ble with the observation that several of the binaries we have 
discovered have short orbital periods corresponding to val- 
ues of log(a/au) ^ —0.5. If it were the case that the bina- 
ries we have detected are biased towards \og{a/au) ^ —0.5 
then the values of Pcmpiricai we calculated for stars in the 
faint sample would have been very small. In this case, the 
probability distribution for the binary fraction calculated in 
section IT2l would have been consistent with the hypothesis 
/faint = /bright - In fact, this hypothcsis can be rejected with 
90 percent confidence, so the binaries we have detected are 
not biased towards log(a/au) ^ —0.5. 

The lack of binaries in our sample with masses ^ O.IMq 
is consistent with the overall picture established from exist- 
ing surveys that the binary fraction for VLMS and BDs 
acros s the full range of orbital separat ions is 20-25 per- 
cent (|Allenll2007l : iBasri fc Reinerslbood ). However, the bi- 
nary properties for VLMS and BDs is poorly characterized 
in the range a ^ 2 au that is below the detection limits 
of high angular resolution imaging. Some VLMS and BDs 
show apparent RV shifts comparable to the orbital speeds 
expected at these separations, and a complete spectroscopic 
orbits has now been publi shed for one such star (Cha Ha 8, 
IJoergens fc Muelleij l2007h . It is still possible given all the 
observations to-date that there is a large population of 
VLMS/BD binaries with a « 1 an. If such a population ex- 
ists then it may be possible to reconcile the binary fraction 
in existing surv eys with the h i gh bin ary fraction inferred for 
the Pleiades bv lLodieu et all l|2007h . 

The spectra we have used for this study were obtained 
with the aim of measuring the binary properties of the low 
mass stars and brown dwarfs in the A Ori and a Ori clusters 
but there is clearly useful information to be obtained about 
the properties of the clusters themselves from the spectra we 
have obtained. For example, the distribution of EW(Nal) 
shown in Fig. [T^] clearly shows that the stars in A Ori have 
higher surface gravities on average than stars in Group 2 
of the a Ori cluster. This shows that the A Ori cluster is 
older than the Group 2 of the <j Ori cluster and, perhaps, 
has a similar age to Group 1 of the a Ori cluster. This is in 
general agreement with the ages we have adopted for these 
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populations. A more detailed analysis is beyond the scope 
of this paper. 



6 CONCLUSIONS 

We have conducted a large radial velocity survey for spectro- 
scopic binary stars among low mass stars and brown dwarfs 
in the young clusters around A Ori and a Ori. Wc have iden- 
tified 196 members of these clusters based on their radial 
velocity and spectral type. Of these, 6 are SB2 binaries 
and 6 are SBl binaries. All the spectroscopic binaries we 
have detected arc brighter than Ic = 16.6 (mass ~ 0.12Mq). 
We conclude that the frequency of spectroscopic binaries in 
these clusters among very low mass stars (mass < O.IMq) 
and brown dwarfs is significantly lower (< 7.5 percent) than 
that for more massive stars (9 ± 2 percent) . The change in 
binary properties with mass that we have discovered may 
be due to a change in the total binary frequency with mass 
or a change in the period distribution of the binaries with 
mass or both. 

The number of SB2 binaries in this sample suggests 
there may be a peak in the mass ratio distribution for spec- 
troscopic binaries in these clusters near g « 1. There is also 
clear evidence from the properties of the SBl binaries that 
the mass ratio distribution for spectroscopic binaries in these 
clusters is broad, extending down to g w 0.25. 
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